KOBE-FHD-99-05 
December 1999 



Two— Spin Asymmetry for i/j f Photoproduction 
with Color— Octet Mechanism 



K. Sudohfl and T. MoriiQ 

Faculty of Human Development and Graduate School of Science and Technology, 
Kobe University, Nada, kobe, 657-8501, Japan 



Abstract 

We studied the photoproduction of ip' in the forward regions in polarized 7p 
collisions at relevant HERA energies. We found that this reaction is very effective to 
test the color-octet mechanism which is based on the NRQCD factorization formalism. 
Furthermore we found that the value of the NRQCD matrix elements can be severely 
constrained by measuring the two-spin asymmetry, though the process depends on the 
polarized gluon distribution Ag(x). 
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Heavy quarkonium productions and decays have been traditionaly calculated in 
the framework of the color-singlet model |]J. However, it has been reported that the 
color-singlet model cannot explain the experimental data of J /if) and if)' hadropro- 
ductions 0; the cross sections of prompt J /if) and if)' production in unpolarized pp 
collisions predicted by the color-singlet model were smaller than the Tevatron data by 
more than one order of magnitude ||. Furthermore, in the photoproduction of J/if>, 
the color-singlet model cannot satisfactorily reproduce cross sections for inclusive J/ if) 
production at the recent HERA experiment M even including next-to-leading order 
corrections ||. In order to solve these serious problems, a new color-octet model has 
been advocated by several people f| as one of the most promising candidates that 
can possibly remove such big discrepancies between the experimental data and the 
theoretical prediction by the color-singlet model. 

A rigorous formulation of the color-octet model has been introduced based on 
an effective field theory called nonrelativistic QCD (NRQCD), in which the 0(v) cor- 
rections of a relative velocity between the bound heavy quarks can be systematically 
calculated 0. The NRQCD factorization approach separates the effects of short dis- 
tances that are comparable to or smaller than the inverse of heavy quark mass, from 
the effects of longer distance scales of hadronization. A heavy QQ pair is first pro- 
duced in a virtual color-octet intermediate state of the NRQCD higher Fock state, 
and then hadronizes into a detected color-singlet particle via the emission or absorp- 
tion of dynamical gluons. Production cross sections of heavy quarkonium H can be 
factorized into a product of a short distance coefficient C n which can be computed 
using perturbative QCD, and a long distance part (O^ ( 2S+1 Lj)) which is described 
by nonperturbative NRQCD matrix elements whose values should be determined from 
experiments or lattice gauge theory, 

£ da(AB ^H + X) = ±dTY: C n (0»C s+1 Lj)), (1) 

X * n 

where $ and V show a flux and a phase space factor, respectively. The label n denotes 
color and angular momentum numbers. In this factorization approach, there are two 
long distance parameters which are essentially the vacuum expectation values of the 
color-singlet and -octet NRQCD matrix elements (0^ 8 ), whose relative importance is 
determined by the velocity scaling rules j5|. Physics of the color-octet model is now 
one of the most interesting topics on heavy quarkonium productions at high energy. 
However, although the color-octet model is quite succesful in explaining the Tevatron 
data, it looks problematic for the process 7 + p — > J /if) + X; using the long distance 
matrix elements extracted from the Tevatron data, the color-octet model largely over- 
estimates the HERA data 0] . Since the color-singlet matrix elements is related to the 
radial wave functions at the origin, their values are calculable using potential models 
|T0|j . On the other hand, the color-octet matrix elements can be extracted only from 



experiments at present. However, the present uncertainties on the color-octet matrix 
elements are still large and the discussion seems controversial. In order to confirm 
the validity of the NRQCD factorization approach, the universality of long distance 
parameters must be established for different processes with acceptable experimental er- 
ror. So far, several processes have been already suggested for testing this factorization 
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Figure 1: Feynman diagrams for the leading order of the color-octet subprocess 7+p — > 
(cc) octe £ and of the color-singlet subprocess 7 + p — > (cc) sing i et + g. Initial 7 and proton 
are polarized. 



approach, such as polar angle distributions of J/ip in e + e annihilation into J /if) + X 
|TT| , Z° decays at LEP |Q, leptoproduction of J/^ ]13[, and so on. 

In this paper, as another test of the color-octet model, we propose the ip' pho- 
toproduction at small-p^ regions, 



7 + p — > i>' + x, 



(2) 



which might be observed in the forthcoming polarized HERA experiment, where the 
arrow attached to particles means that these particles are polarized in a parallel or 
antiparallel direction to the running direction of each particle. In this work, we do 
not consider the contribution of elastic diffractive mechanisms and other reducible 
background processes, since they can be eliminated with appropriate cuts |14j. If the 
color-octet mechanism works, the process is dominated by the color-octet subprocess 
7 + g — > (cc) oc t e t at small-pr regions. In addition, there can be the contribution from 
the conventional color-singlet subprocess 7 + g — > (cc) sing i e t + g which appears only as a 
second order of the strong coupling constant a s . Feynman diagrams for these subpro- 
cesses are illustrated in Fig. [I]. The main contribution to the leading color-octet mech- 
anism comes from the color-octet ^q 8 ^ and 3 Pj 8 2 02 states, whereas the color-singlet 

mechanism originates from 3 S[^ state. The process involves the NRQCD long distance 
parameters corresponding to such states. In particular, for the color-octet machanism 
the cross section is propotional only to the linear combination of the NRQCD long 
distance parameters, {0% ( 1 5'o)) and (Of ( 3 -Po)), owing to the NRQCD spin symmetry 
relation, as described later. We show that to measure the spin-dependent cross section 
and two-spin asymmetry for this process is very effective not only for testing of the 
color-octet model but also for extracting the values of the color-octet long distance 
matrix elements. Furthermore, since the process is dominantly produced in photon- 
gluon fusion, the cross section must be sensitive to the gluon density in a proton and 
thus one can get good information on the spin-dependent gluon distribution function 
by analyzing this polarized process. A related subject has been recently investigated 
by Japaridze et al. |TJ|]. They have calculated the two-spin asymmetry of J/ip photo- 
production at large-pr regions in polarized 7^ scattering and discussed the sensitivity 
of the long distance parameters to the asymmetry. They have insisted that the process 
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Figure 2: JTie spin-independent total cross section as a function of ^/s. The solid line 
represents the sum of the color-sinlet and octet contributions, while the dashed line 
represents the color-singlet contribution only. 



is effective for testing the color-octet model. 

The spin-independent total cross section due to these mechanisms for small-pr 
regions is given by 

o-(lP -> ij/ + X) = a^p -> ip' + X) + o"i(7p -»• ^' + X) 



2S+1 r (8) 



7r 2 a e me 2 



+ 



dt 



)(Of( 2S+1 Lj)) 
+d-( 1 g^cc[ 3 S[ 1) ]g)(Of'(*S 1 ))) 

dxg(x) (na s 5(Aml — s)6 
64a^|%(0)| 2 s 2 (s - 4m 2 ) 2 + i 2 (i - Ami) 2 + u 2 (u - Am, 



2\2 S 



3s 2 



Am 2 ) 2 (t — 4m 2 ) 2 (w — 4m 2 )^ 



(3) 



where g(x) is the unpolarized gluon distribution function, and s, t and u are the usual 
Mandelstam variables for the subprocess. The labels 1 and 8 for the cross sections 
denote the contribution from the color-singlet and -octet state, respectively. The sum 
in the first term is taken over 3 Pq 2 an d 1 Sq states and is the linear combination of 
color-octet matrix elements for these states, 
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mi 



(4) 



R^i(0) is a radial wave function at the origin and is related to the color-singlet matrix 
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Figure 3: The spin-dependent total cross section with the parameter 0/0 = 3.6 as a 
function of y/s. The solid, dashed and dotted lines show the case of set A of GS96 fiH[] , 
set B of GS96 ftF% and the 'standard scenario' of GRSV96 jlffi , respectively. Upper 
bold lines represent the sum of the color-singlet and -octet contributions, while lower 
lines represent the color-singlet contribution only. 



element as 

(0f( 3 5i))~^|^(O)| 2 (l + ^ 2 )), (5) 

which can be well determined from the leptonic decay width of r/>', though its value 
cannot be calculated perturbatively. 

The spin-dependent cross section can be obtained by replacing the unpolarized 
subprocess cross section by the polarized one, and furthermore by the following re- 
placement, 

g(x) — > Ag(x) (A means "polarized".), (6) 

e - e = (of'CSo)) - — 2 {of ( 3 p )>. (7) 



In the present calculation, we used the GS96(set A and set B) fl7j and GRSV96 
IE| parametrizations for the polarized gluon distribution and the GRV95 parametriza- 



tion ]19| for the unpolarized one. Color-octet matrix elements were taken from the re- 



cent analysis on charmonium hadroproduction data: (0% ('So)) + (7/ml){Og ( 3 -Po)) 
5.2 x 10- 3 [GeV 3 ] and \{Of ( X S Q )) + (l/m 2 c )(Of ( 3 P )> « (5.9 ± 1.9) x l(T 3 [GeV 3 ] 



|21j|, which lead to 

= (Of ('So)) -^(Ot'( 3 Po)) 
O" (Of ('So)) + ^ (Of ("P )) 



3.6 ~ 8.0. 
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Calculated cross sections are presented in Fig. |2| and Fig. We see that the color-octet 
contribution is larger than the color-singlet one by one and two order of magnitude 
for the unpolarized and polarized cross sections, respectively. It is remarkable that the 
difference due to different gluon polarization models is not so large. 

Now we move to the analysis on a two-spin asymmetry for ip' production in the 
polarized reaction defined by 
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LL 



[da ++ — da.\ + da — do 

[d<j ++ + da + - + da + dcr_ H 

dAa dAa$ + dAa\ 
da da 8 + da\ 



(9) 



where da + ^, for instance, denotes that the helicity of photon is positive and the one of 
proton is negative. From Eq. |J we see that if only the color-octet contribution which 
is the lowest order process in a s works, the asymmetry An can be written by a simple 
formula 

^ dAa 8 Ag(x,Q 2 ) 

das g{x,Q) © 

which is just a product of the ratio of polarized and unpolarized gluon densities and 
the one of color-octet matrix elements. 

Taking Q 2 = Am 2 with a charm mass m c = 1.5GeV, the calculated An at 
relevant HERA energies are presented in Fig. |4j. As shown in Fig. ||, if the color- 
octet process works, then the An becomes quite large in the rather smaller yfs region, 
comparing with the one for the color-singlet mechanism only. The difference of An 
due to the color-octet and -singlet mechanism is larger than the uncertainties due to 
the polarized gluon distribution functions. Hence we can sufficiently test the color- 
octet contribution in this reaction. Furthermore, since the An strongly depends on 
the value of 0/0, one can constrain its magnitude from the value of An as follows; if 
we take the GS96 or GRSV96 parametrization which are widely used, the maximum 
value of Ag(x)/g(x) becomes roughly 0.35 for GS96 and 0.2 for GRSV96 and then, 
with this value on the ratio of gluon distributions, we can constrain the maximum value 
of the ratio of NRQCD matrix elements as 

| < 5.0 forGS96, (11) 
< 2.9 for GRSV96, (12) 

from the requirement that the An should be less than 1. Actually, the uncertainty 
of matrix elements seems to be larger than that of gluon distribution, because the 
value of matrix elements obtained from the Tevatron data which we used here does 



not include the contributions of higher order QCD corrections. In Ref. ||22|| , Kniehl 
and Kramer have approximately taken into account the effect from higher order QCD 
corrections due to multiple-gluon initial state radiation and improved the values of long 
distance parameters for J/ip meson as smaller. They have insisted that it is possible 
to explain both Tevatron and HERA data by using such a 'small' set of long distance 
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Figure 4: The two-spin asymmetry All{iP ip'X) with the parameter 0/0 = 3.6 as 
a function of ^/s. Various lines show the same as in Fig. [|. Upper bold lines represent 
the color-singlet plus octet contribution, while lower lines represent the color-singlet 
contributions only. 



parameters. In any case we can say that our process is very effective not only for 
testing the NRQCD factorization approach but also for constraining the value of long 
distance parameters, though the result depends on the polarized gluon distribution 
function. On the contrary, if the NRQCD factorization approach is confirmed enough 
with the long distance matrix elements with acceptable theoretical and experimental 
uncertainties, we can get good information on the polarized gluon distribution function 
in rather smaller y/s regions. 

Finally, let us discuss the sensitivity on our results. In order to examine the 
experimental feasibility of the forthcoming HERA experiments, we have estimated 
the experimental sensitivity of the All for 100-day experiments at various y/s in 



the manner of Nowak ||23|| , using the expected data of beam or target polarization 
(Pb, Pt ~ 70%), the integrated luminosity (C • T ~ 66pb _1 ), and the combined trigger 
and reconstruction efficiency (C ~ 50%) together with the value of unpolarized total 
cross sections. As a result we found that the experimental sensitivity SAll is order 
of magnitude ~ 10~ 3 , which is very small. Hence our predictions are expected to be 
actually tested in the future polarized HERA experiments. 

In summary, to test the color-octet model we have proposed the photoproduction 
of if)' at small-pr regions in polarized 7p scattering which might be available in the 
forthcoming polarized HERA experiments. We have calculated two-spin asymmetry 
All for various parameter regions O/O = 3.6 ~ 8.0, and found that the All becomes 
quite large in the regions \/s = 10 ~ 20 GeV. Therefore we can sufficiently test the 
color-octet model in this process. In addition, the measurement of All is very effective 
to severely constrain the value of NRQCD matrix elements, though it depends on the 
polarized gluon distribution Ag(x). 
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